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Simple field-theoretical approach to critical phenomena is described. In contrast to the Wilson's 
theory, a description in real 3-dimensions space is used. At the same time the described approach is 
not the same as Parisi's. Used subtraction scheme is different from the one used by Parisi, but the 
main point is that we treat changes of T — Tc as explicit perturbations. By such an approach not 
only the critical domain but also the crossover from critical domain to domain of Landay theory can 
be described very simply. Consideration is restricted to scalar 4 model in one-loop approximation. 
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The problem of critical phenomena description is known for a long time and was successfully solved in 70's mainly 
by Wilson^. Wilson's theory is field-theoretical one and essentially use formal expansion in series in e, so called 
O ■ e-expansion. Here e is deviation of space dimensions from a case of 4 dimensions. 

Besides Wilson's approach using e-expansion yet another field-theoretical approach working directly in the real 3- 
dimensional space is also known. This approach proposed by Parish has been used in calculations of critical exponents 
+Jj ' of the scalar model up to five-loop approximation^. Similar approach was proposed by Ginzburg^. 

Up to now both these approaches have been developed in detail and now they form the standard basis of the critical 
r/3 ■ phenomena theory (see£~— and references therein). 

However yet another field-theoretical approach to critical phenomena is possible. This approach is described in this 
paper. As in Parisi's approach we use renormalization group (RG) in real 3-dimensional space. But in contrast to 
latter we use point of renormalization at arbitrary momentum. Such an arbitrary point lets us to describe a massless 
field as well. Besides, we treat the changes of T — Tc as explicit perturbations by additional term in the Lagrangian. 
Although such details are not critical for physical meaning of a theory, these changes allow to simplify essentially 
the physical sense of a theory. Besides that, such an approach allows to describe very simply not only the critical 
domain but also the crossover from critical domain to the domain of Landay theory. It is worth to point out that the 
aim of this paper is not to calculate critical exponents with high accuracy but to draw simple and intuitively clear 
physical picture of a critical phenomena. This is why the consideration is restricted to scalar tp 4 model in one-loop 
' approximation and consideration is as simple as possible. 
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l/-) ; II. PERTURBATION THEORY AND DIAGRAMS 

In this section we describe very briefly perturbation theory to introduce the notations. The details can be found in 
1 any quantum fields textbook, say by Ramond^. 

Statistical physics of thermally fluctuating scalar field 4> with interactions described by term 4 is defined by 
. . ■ partitional function which can be represented as path integral 



Z{ J) = J e~ s ^ J) V<i> , (1) 
3(4>, J) = S (<f>, J) + Sj(^) = J £ d 3 x + J £/d 3 x , (2) 

^o = i(V0) 2 + ^ 2 -J0, C I = ^4> i . (3) 

Generally speaking it should be additional T _1 in the Gibbs exponent in ([I} and some coefficient / in (V</>) 2 /2 in 
©. But one can redefine field (j) — > T -1 / 2 / 1 / 2 ^ so that these coefficients disappear. Given such field redefinition To , 
go and J are also redefined: 

r -> ro/- 1 , (4) 
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where "action" is 



and "Lagrangian" is 



9o 9oTf~ 
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(5) 



J^JT- 1 ' 2 /- 1 ' 2 . (6) 

It is easy to see that variational derivatives of Z(J) by source field J(x) generates non-normalized correlators 
G = Z(0)((f>(xi) . . . 0(x n )) which can be named Green's functions in analogy with quantum field theory. To make 
path integral calculable, the exponent with the interaction term is expanded in series: 



Such a representation yields the well known perturbation theory which can be represented by Feynman diagrams. 
These diagrams consist of solid lines and elementary vertexes. A solid line corresponds to the propagator (two-point 
Green's function) with Fourier transform (k 2 + To) -1 . Elementary vertex is — go/(4\) . 

Further one can introduce connected Green's functions Gc and so called one participle irreducible Green's functions 
r in a common way. Self-energy function (mass operator) X(fc) in a diagrammatic form is presented as a sum of two- 
point diagram (without external tails) which can not be divided in two parts by removing of one line. With self-energy 
function one can easily get the exact perturbated propagator: 

Other details see in the textbook cited above. 



III. RENORMALIZATION OF INTERACTION PARAMETER 



Let us describe 4-point one-participle irreducible Green's function . It is defined by diagram: 




Analytical representation of this diagram is: 

-r 4 (pi,p 2 ,p3,P4) = .9o - gl I i < i) > ( 9 ) 

where /(q) are loop integrals different from each other by permutations of diagram's tails; q is the momentum which 
flows throw a loop, it can be found from tail momenta p a by momentum conservation law. 

After that we renormalize the perturbation theory. Usually renormalization is used when loop integrals diverge. 
In three dimensions </> 4 is super-renormalizable theory, there is the only divergent diagram and it is different from 
the shown above. But anyway it is possible to renormalize the convergent diagrams also. Such a renormalization 
improves the perturbation theory. 

To renormalize perturbation theory we divide parameters into couples of terms tq = r + r ct , go = g + g ct and 
treat the term of Lagrangian with r c t as perturbation. Besides that, 4 -terms with g and g c t are treated as separate 
perturbations. Because the new perturbation series is series in g , counterterms r ct ~ g and g c t ~ g 2 at least. Generally 
speaking, the gradient term (V^) 2 /2 should be divided into a main term and a counterterm as well. But it is not 
necessary in one- loop approximation of 4 -theory and we do not describe it. 

According to the above, renormalized equation © up to g 2 is: 

- r 4 (pi,p 2 ,P3,P4) = g- g 2 ^2l(q) +g c t- (10) 

Certainly there is a continuum of variants to divide bare interaction constant go into renormalized interaction param- 
eter g and counterterm g ct . Thus we should add some renormalization conditions to choose one of them. 

Without perturbation theory it does not matter what variant to choose. But if one uses perturbations then some 
variants are much better then the other. Obviously perturbation expansion works well when third term in r.h.s. 
of (jlOp approximately compensates the second term. Thus it is convenient to choose renormalization condition as 
a condition of full compensation of the second and the third term when the tail momenta are equal to the some 
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four vectors, so called renormalization point vectors. Then the perturbation expansion will work well when the tail 
momenta are close to renormalization point vectors. 

Now we should choose the normalization point. To parametrize four vectors by one parameter we introduce for 
unity vectors n a and define normalization point as p a = /in a where fj, is renormalization parameter. It is convenient 
to choose vectors n being inscribed in cube symmetrically as shown on FIGQ] Obviously with such vectors n a vectors 
p a = pn a obey momentum conservation low. 




FIG. 1: Symmetric set of vectors n a . 



Further we restrict the consideration of r^pi, P2, P3, P4) to a case when p a = pn a . Such a T4 will be denoted as 
^i(p) ■ For such a case q is the same for all permutations of tails and we can write: 



- r 4 (p) = 1 

With elementary geometry one can find q 



3 2 1 
2 9 (2tt) 3 



d 3 k 



(fc 2 +r)[(k + q) 2 +r] +5ct - (U) 

(2p) I \f?> . According to said above we define renormalized g by condition: 

-T 4 (p)\ p ^=g. (12) 

Loop integrals similar to integral in (jlll) are usually evaluated by a well known Feynman trick. Specifically this 
integral also can be calculated by elementary methods in cylindrical coordinates with z-axis along vector q . Anyway 
integration and fixing of g ct by (fL2"| yields: 



r<i(p) =9- ff 



3V^ 

167T 



• arctan 



/3r 



arctan 
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(13) 



Actually equations (|12[) and (|13[) define how renormalized interaction parameter g changes when the scale of 
renormalization point /i changes. One need only substitute the new scale \j! as p in (113[) . But it is worth to remember 
that (fT3"|) works well only when p is close to fi . Thus it is much better to describe large change of ijl as step-by-step 
small changes of fi . If each step is infinitesimal then this yields a differential equation. This is the key idea of 
renormalization group method. Substituting p = fi + dfi in (113[) with simple algebra we get: 



dgy 
d/j 
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3\/3 J_ 

167T /i 2 



arctan 



/' 



M/v3r 



1 + (/VV3r) 5 



(14) 



where the subscript fi stresses that g is a function of fi . This equation is one of the key equations of our description. 
When t = it becomes more simple: 



dg^ _ 2 3\/3 
d/j ~ 9>1 32 



1 



(15) 



IV. RENORMALIZATION OF "MASS 



The parameter r corresponds to the square of the participle mass in quantum field theory. This parameter also 
undergoes renormalization. In one-loop approximation this renormalization represented by single diagram for self- 
energy function: 
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It is essential that in one-loop approximation E is momentum independent. In our theory r is /i-independent constant 
while one-loop approximation is used. 

It is worth to note that in the common theories r is /i-dependent nevertheless. This is reached by very refined trick 
which we do not use. Instead, we use another more physically clear trick which is described below. 



V. "RUNNING" INTERACTION PARAMETER 



In this section we solve equation (fT4"]) but before that it is useful to describe solution of (TT51) . Obviously, when 
/i >> s/t the solution of (TT4|) is close to the solution of (|15l) . while the latter is much more simple. Solution of (fTS"|) 
is: 

8= aTcH> (16) 

where A = 3V3/32 w 0.1624 , and C is an integration constant. One can see easily that if /i — > oo then g — > C _1 . 
Thus at large momentum g is constant. Let us denote this constant as g^ . Then 

Otherwise if /i — > then 

5 ( M ^0) = ^«6.158 M . (18) 

Equation (1181) shows a very important fact: the behaviour of g in IR limit is universal, it does not depend on g^ . 
This means that if <fi describes some field in real system, say in crystal, then interactions of its long-wave thermal 
fluctuations does not depend on interactions at atomic scale. Now we see this fact for r = (i.e. for a system at just 
critical point), further we will see this for g^ »t^0 as well. 

Exact solution of (fT4|) yields not too complex result: 

9 = 7= 7= . (19) 

[(3V3)/(16tt)] arctan(^/V3r) + C/i 

where C is an integration constant. We see again that g — > const = C _1 = g^ while ^ — > oo . So we can rewrite the 
equation in the form: 

g(p) = ?=-^ /=- • (20) 

1 + .9ooM 1 [(3v3)/(167r)] arctan(^/V3r) 

Some essential point should be discussed here. Describing solution of differential equation mathematically, one 
may think that the integration constant C and the parameter g^ can depend on r . But by physical reasons this 
parameters should be r-indepedent. Indeed, one can see from (|13l) that if fj, — > oo, then g c t — > . Thus g^ actually 
is bar parameter go which obviously does not depends on r . 

Further, if fi = then: 

.9(^ = 0) = -——^——^. (21) 
1 + 3oo[3/(167t)Jt V* 

One can see very easily that if r << g^ then the behaviour is universal again: 

.9(^ = 0) = ^V?« 16.755V?. (22) 



Besides if g^ » fx >> y/r then 



32 

9 (/x) = ^ M « 6.158/x. (23) 



This equation is also universal and it is the same that was derived above for the massless (r = 0) field. 
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VI. PERTURBATION OF PROXIMITY TO CRITICAL POINT 



Experimentally, the proximity of a system to the critical point is usually controlled by temperature. Thus a direct 
way to describe system reaction to the changes of proximity to critical point requires to find a reaction to g Q changes 
(see equation ([5]) ). But if T — T c « T c then it does not matter what to change, T or T c ■ This is why the proximity 
to critical point can be controlled by tq changes, this is more convenient technically. 

We describe T — Tq changes as an explicit perturbation. Thus we add a perturbation term tXo4> 2 /2 to the 
" Lagrangian" . Here the parameter t controls the changes of T — Tc . Why a parameter Ao is added will be clear 
further. 

In reality, t is spatially homogeneous. However, perturbation with spatially homogeneous t can not be renormalizcd 
directly. This is why we generalize the theory to spatially inhomogeneous i(x) which will be changed to homogeneous 
t at the end of calculations. Actually we add yet another source held interacting with <f> non-linearly. 

Perturbation i(x)Ao</> 2 /2 generates additional 3-tail vertex (— Ao/2 in the analytical representation). One of its 
tails is £(x) which will be depicted in the diagrams as a waving line. As any vertex, A- vertex can be renormalized. To 
renormalize it we divide it into renormalized vertex A and countcrterm X ct where X c t ~ g at least. Further descriptions 
are generally the same as the renormalization of g . One-loop correction to T 3 is represented by diagram: 




+ 



+ 



which analytically yields: 




(24) 



where p is waving tail momentum. Renormalization conditions we choose similar to condition for g : 



r 3 (p)| 



= A. 



(25) 



Integral in (|24p is just the same that in (fTTj) . So with full analogy to above descriptions we derive: 




(26) 




(27) 



For t = or /j, » \fr the last equation can be simplified: 



d f i M 16 ' fi 2 



(28) 



VII. 



"RUNNING' 



A 



A(/i) can be found easily when r = or // >> y/r . Using (ITTl) and ([28|) we obtain: 




(29) 
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where a = 2/(3^3) « 0.385, C is an integration constant. Obviously if —> oo then A — > C — const . Thus the 
equation can be rewritten as: 

Here the parameter X^, is r -independent by the same reason as . When fi « g^ the equation can be simplified: 

x(p)= [ 9oo (31) 

If fi ~ •y/r" or fi < y/r then the differential equation (1271) apparently can be solved by numerical methods only. But 
at least in critical domain it is easy to get universal correlation between A(/x = 0) and A(/Lt*) , where //* ~ ^/r, say 
//* = 10-y/r . Further A(/z*) can be evaluated by (|30|) or ([3T]) . 

If < [i* = lOy/r and fi* << g^ then (|20|) can by simplified as 

ffto = P ^ P= ■ (32) 
3v 3 arctan(/z/v 3r) 

With such g(/j,) equation (f27|) yields: 
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dA = ^ > [arctan(x/2) - (s/2)/[l + (x/2) 2 ]] ^ 
A 7 x arctan(a;/v3) 

A(0) 

where substitution x = \il\fr used. 

In (1331) r.h.s. is a universal constant which can be evaluated numerically. For our purpose however it does not 
matter what its value is. The only fact we need is that according to ([55]) 

A(0)~A(/i*). (34) 

With (|3ip this correlation yields in critical domain: 

A(0) ~ W 2 . (35) 



VIII. CRITICAL EXPONENTS 



The last equation of the previous section allows us to derive critical exponents very simply. If the system is off the 
critical point by some non-zero r , then, in order to bring it to the critical point, the perturbation of the self-energy 
function E(fc = 0) should be equal r. On the other hand, in the first order of perturbation theory £ = — A(0)AT, 
where AT is the temperature shift (negative to shift down to critical temperature). Thus we have a self-consistent 
equation: 

r = -A(0)AT. (36) 

Using (|35l) we instantly obtain: 

T ~T a/2 (T-T c ). (37) 

An obvious transformation yields: 

r ~ (T - T c y , (38) 

where 



2 

1= —a 



(39) 
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There is no field renormalization in one-loop approximation. This is why the Fisher exponent ry is zero and the 
propagator is G(k) = (k 2 + r) _1 . Thus the critical exponent for correlation length is 

"=i = irb- (40) 

Numerical estimations yields values 7 = 1.24 , v = 0.62 which are close to the known precise values (see£) . 

Derivations in this section above may look doubtful because A(0) changes while t changes. Using of the first order 
perturbation may also seem doubtful. But we can justify the result by consideration of infinitesimal change of r with 
infinitesimal change of temperature. Such an approach yields: 

ch = \(Q)dT ~ T a/2 dT . (41) 

Integration of this differential equation yields exactly the same result that we have obtained before. 



IX. CROSSOVER TO LANDAU THEORY 



The approach considered here allows us to determine not only critical behaviour but also crossover to Landau theory 
behaviour when r is large. To describe this crossover one should find A(0) as a function of r without approximations. 
Using (HOI) and {H} we obtain 



A(0) = A M e~ f M 



where function F(y) is defined as 



00 

, f arctan(a;/2) - (x/2)/[l + (x/2) 2 } dx 

F (y> = y / ; — ~, — : — ' 

J x + ya~ L arctan(a;/ yo) 



(42) 



(43) 



V = 9oo/ (87T-v/t) , x = [J,/ \fr . A dependence of A(0)/Aoo on y = goo/ (&tt*/t) can be calculated by numerical methods. 
The result is shown on FIGJ2J 



A(0) 




0.0001 



FIG. 2: Dependence of A(0)/Aoo on y — g 00 /(8TT^/r) calculated numerically. 



From FIG(2]one can see that if \[t << goo/ (Sir) i.e y » 1 then A(0) depends on r as described in previous sections. 
This behaviour corresponds to critical domain. Otherwise if y/r » g^f^n) i.e. y « 1 then A(0) ps Aqo = const. 
This means that in this domain dr ~ dT . Such a behaviour corresponds to Landau theory. One can see also that 
condition of crossover r ~ 5oo/(87r) i.e. y ~ 1 coincides, at least up to some coefficient, with well known Levanyuk 
criterion. 
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X. CONCLUSION 



In this paper a very simple approach to critical behaviour description has been presented. Generally this approach 
is close to the common field-theoretical theories using renormalization group method. Nevertheless the presented 
approach has at least one advantage: it allows one to draw a simple, intuitively clear physical picture of a critical 
phenomena. 

Physical reason why non-trivial exponents appear is that the system's sensitivity to T — Tq variations, described 
by A- vertex at zero momentum, changes as system goes toward to critical point. If one takes in account this fact, 
then reasonable values of non-trivial exponents can be obtained by very simple calculations. Moreover, far from 
critical point A- vertex ceases to depend on correlation length (which is ~ l/^/r). In this way, far from critical point 
the behaviour of a system becomes identical to the one described by Landau theory. Such a crossover can also be 
described in a simple way. 
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